Spatially Resolved Raman Spectroscopy of Single- and Few-Layer Graphene 
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We present Raman spectroscopy measurements on single- and few-layer graphene flakes. Using 
a scanning confocal approach we collect spectral data with spatial resolution, which allows us to 
directly compare Raman images with scanning force micrographs. Single-layer graphene can be 
distinguished from double- and few-layer by the width of the D' line: the single peak for single-layer 
graphene splits into different peaks for the double-layer. These findings are explained using the 
double-resonant Raman model based on ab-initio calculations of the electronic structure and of the 
phonon dispersion. We investigate the D line intensity and find no defects within the flake. A finite 
D line response originating from the edges can be attributed either to defects or to the breakdown 
of translational symmetry. 

PACS numbers: Valid PACS appear here 
Keywords: Raman mapping 
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The interest in graphite has been revived in the last 
two decades with the advent of fullerenesi and carbon 
nanotubes. 2 However, only recently single- and few-layer 
graphene could be transferred to a substrate^ Transport 
measurements revealed a highly-tunable two-dimensional 
electron/hole gas of relativistic Dirac Fermions embed- 
ded in a solid-state environment 4£ Going to few-layer 
graphene, however, disturbs this unique system in such 
a way that the usual parabolic energy dispersion is re- 
covered. The large structural anisotropy makes few-layer 
graphene therefore a promising candidate to study the 
rich physics at the crossover from bulk to purely two- 
dimensional systems. Turning on the weak interlayer 
coupling while stacking a second layer onto a graphene 
sheet leads to a branching of the electronic bands and 
the phonon dispersion at the K point. Double-resonant 
Raman scattering* which depends on electronic and vi- 
brational properties turns out to be an ingenious tool to 
probe the lifting of that specific degeneracy. 

We report on Raman mapping of single- and few-layer 
graphene flakes resting on a silicon oxide substrate. Lat- 
eral resolution of 400 nm allows to address neighboring 
sections with various layers of graphene down to a sin- 
gle graphene sheet, previously determined with the scan- 
ning force microscope (SFM). We find that the integrated 
G line signal is directly correlated with the thickness of 
the graphitic flake and is shifted upward in frequency 
for double- and single-layer graphene compared to bulk 
graphite. The mapping of the peak width of the D' line 
shows a strong contrast between single- and few-layer 
graphene. Such a pronounced sensitivity to the tran- 
sition to the very last layer offers an optical and non- 
destructive method to unambiguously detect single-layer 



graphene. In addition, we locally resolve the structural 
quality of the flake by investigating the D band, which 
is related to elastic backscattering. The map of the in- 
tegrated D line signal of a graphitic flake with double- 
and single-layer sections shows that the inner part of 
the flake is quasi defect free, whereas edges and steps 
serves as scatterers. Finally, we explain the splitting of 
the D' line as a function of the number of graphene layers 
within the double-resonant Raman models The compar- 
ison between experimental data and theory confirms the 
qualitative validity of the double-resonant Raman model, 
but demonstrates quantitative differences between theory 
and experiment. In particular, the model, when based 
on first-principles calculations, predicts a much smaller 
splitting of the peaks. 

The graphite films were prepared by mechanical ex- 
foliation of highly oriented pyrolytic graphite (HOPG) 
and subsequent transfer to a highly doped Si waver with 
300 nm Si02 (atomic oxidation process) cap layer The 
combination of optical microscopy using phase contrast 
and SFM makes it possible to locate flakes with various 
thicknesses down to a monolayer with lateral extensions 
in the micrometer range. The Raman spectra were ac- 
quired using a laser excitation of 532 nm (2.33 eV) deliv- 
ered through a single-mode optical fiber, whose spot size 
is limited by diffraction. Using a long working distance 
focusing lens with a numerical aperture NA=0.80 we ob- 
tain a spot size of about 400 nm. With a very low incident 
power of 4-7 fiW heating effects can be neglected. 

The Raman spectrum of graphite has four prominent 
peaks (Fig. |3- for a recent review see Ref. |8|). The peak 
around 1582 cm -1 , commonly called G line, is caused by 
the Raman active Ei g phonon, (in-plane optical mode) 
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FIG. 1: (a),(b) SFM micrograph and cross-sectional plot (in- 
dicated with the white dashed arrow; lateral average over 400 
nm) of a few-layer graphene flake with central sections down 
to a single layer. Raman maps (dashed square corresponding 
to the SFM image in (a)) showing (c) the integrated intensity 
of the G line and (e) the FWHM of the D' line. The related 
cross sections (d),(f) are aligned (vertical dashed lines) with 
the height trace. 



close to the T point. The D line around 1350 cm -1 
exhibits two remarkable features: its position shifts to 
higher frequencies with increasing incident laser excita- 
tion energies^ and its relative signal strength (compared 
to the G line) depends strongly on the amount of disorder 
in the graphitic materiali2ii£ The associated overtone D' 
around 2700 cm -1 is pronounced even in the absence of a 
D signal. Finally, the overtone of the G line, the G' line, 
is located at 3248 cm -1 , which is more than twice the 
energy of the G line. The different experimental findings 
related to the dispersive D, D' bands could be explained 
by Thomsen and Reich within the framework of double 
resonant Raman scattering, 6 which was extended to other 
phonon branches by Saito et alM> The electronic and vi- 
brational properties of graphite are dominated by the 
sp 2 -nature of the strong intraplane covalent bonds. The 
relatively weak inter-layer coupling causes the high struc- 
tural anisotropy. Raman spectra for multiple graphene 
layers can be compared qualitatively and quantitatively 
while investigating flakes with sections of various thick- 
nesses. In Fig. ^a) the SFM micrograph of a graphite 
flake with different layers is presented: The bare Si02 (in- 
dicated by '0') is surrounded by single-layer sections with 
steps of up to two, six and four layers. The different step 
heights are clearly depicted in Fig. ^b), where a cross 
section of Fig. ^a) (see white dashed arrow) is shown. 
Scanning the flake and collecting for each spot the com- 




et 0.4 0.8 1.2 1.6 2 1300 1400 



Lateral position (|xm) Raman shift (cm-1) 

FIG. 2: (a) SFM micrograph of a graphitic flake consisting 
of one double- and two single-layer sections (white dashed 
line along the boundaries) , highlighted in the Raman map (b) 
showing the FWHM of the D' line, (c) Raman mapping of the 
integrated intensity of the D line: A strong signal is detected 
along the edge of the flake and at the steps from double- to 
single-layer sections, (d) Raman cross section (white dashed 
arrow in (c)): Staircase behavior of the integrated intensity 
of the G peak (solid line) and pronounced peaks at the steps 
for the integrated intensity of the D line (dashed line), (e) 
Spatially averaged D peak for the crossover from double to 
single layer (•, dashed line) and from single layer to the Si02 
substrate (■, solid line). 



plete Raman spectrum we can subsequently filter specific 
spectral data for spatially resolved data point and con- 
struct false-color 2D plots. In Fig. QJc) the intensity of 
the G peak is integrated from 1537 to 1622 cm -1 . We find 
a remarkable correlation with the SFM graph: Brighter 
regions correspond to thicker sections. The cross section 
in Fig. ^d) shows a step-like behavior, perfectly corre- 
lated with the topographical changes shown in Fig. db) . 
In Fig. □>) we plot the FWHM (full width at half max- 
imum) of the D' line. It shows the narrowing at the 
transition to a single-layer (see e.g. Fig. |3J) and gives an 
evident contrast between single- and few-layer graphene 
sections. The quasi digital change from about 60 to 30 
cm -1 shown in Fig. ^f) suggests that the width of the 
D' line can be used as a detector for single-layer graphene 
resting on a substrate. Raman spectroscopy can therefore 
be used to count the layers of a thin graphite stack and to 
discriminate between single and double layer. Combined 
with the double-resonant Raman scattering mechanism 
an optical setup using light in the visible range turns out 
to be an alternative to scanning force microscopy, which 
requires stacking folds as height references. 

Transport measurements show that the quality of the 
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FIG. 3: Raman spectra of (a) single- and (b) double-layer 
graphene (collected at spots A and B, see Fig. |2Tb11. 
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FIG. 4: (a) Plot of the ratio of the integrated intensities of the 
G and D' peak versus number of stacked layers (average value 
and standard deviation) . (b) G line frequency versus number 
of stacked layers (average value and standard deviation) . (c) 
G peak for HOPG (upper peak), double- (middle peak) and 
single-layer (lower peak) graphene. The vertical dashed line 
indicates the reference value for bulk graphite. 



finite graphitic flakes on the silicon oxide matrix obtained 
with the technique explained above is remarkable: elec- 
tronic mobilities up to 15'000 cm 2 (Vs) -1 were estimated 
from transport experiments 4£ We point out that also the 
Raman spectroscopy reveals quasi defect free graphitic 
sheets via the absence of a D band signal. First exper- 
iments have related the intensity of the D band to the 
structural coherence of the graphite material. In fact it is 
inversely proportional to the crystallite grain size^S The 
appearance of the D band can, however, be related to the 
occurrence of defects and disorder in general, as shown in 
experiments with boron-doped and electrochemically ox- 
idized HOPG>i£ With micro- Raman mapping we are able 
to localize the spatial origin of the defects. From cross- 
correlating the SFM micrograph in Fig. EJa) with the 
Raman map of the integrated D line (1300-1383 cm -1 ) 
intensity in Fig. [^c) we infer directly that the edges of 
the flake and also the borderline between sections of dif- 
ferent height contribute to the D band signal whereas the 
inner parts of the flakes do not. This is somewhat sur- 
prising since for thinner flakes the influence of a nearby 
substrate on the structural quality should be increasingly 
important. In the cross-section Fig. |2{d) we see clearly 
that the D line intensity is maximal at the section bound- 
aries, which can be assigned to translational symmetry 
breaking or to defects. However, we want to emphasize 
that the D line is still one order of magnitude smaller 
than the G line. In Fig. |3{e) spatially averaged D mode 
spectra from the two steps shown in Fig. Efd) are pre- 
sented. The frequency fits well into the linear dispersion 
relation of peak shift and excitation energy found in ear- 
lier experiments. 9 In addition, we find that the peak is 
narrower and down-shifted at the edge of the single-layer 
while it is somewhat broader and displays a shoulder at 
the crossover from the double to the single layer. 

In Fig. 02a) and (b), we compare the Raman spec- 
tra of the double- and single-layer graphene shown in 
Fig. Hfb) and labeled with A and B. The Raman signal 
is significantly altered when peeling off the penultimate 
layer: the G peak decreases strongly in intensity and 
shifts towards higher wave numbers. In connection with 
Fig. ^b) we already stated that the integrated G line 



signal is monotonously increasing with increasing flake 
thickness. In order to compare data of different flakes 
and measurement runs we turn our attention to the ratio 
of the integrated intensities of the G and D' line, plotted 
in Fig. Eta). Most of the changes can be attributed to the 
decrease of the G line, since the spectral weight of the D' 
band changes only slightly. The intensity ratio increases 
almost linearly from one to four layers. In Fig. |l£b) the 
dependence of the G peak position on the layer number 
is investigated. Spectral data of various sections on dif- 
ferent flakes were averaged. The frequency shifts towards 
higher wave numbers at the crossover to the double- and 
especially to the single-layer graphene. However, in the 
case of single-layer graphene, it is accompanied by an im- 
portant statistical spread of the collected data. In Fig. 
^Jc) representative G peak spectra for single-, double- 
layer graphene and HOPG are presented. It is important 
to note that in contrast to the G line, the corresponding 
overtone band, the G' line, does not change its spectral 
position as a function of the number of layers. 

The most prominent difference in the spectra of single- 
layer, few-layer, and bulk graphite lies in the D' line: the 
integrated intensity of the D' line stays almost constant, 
even though it narrows to a single peak at lower wave 
number at the crossover to a single layer (Fig. [3J . The 
width of the D' peak or - at high resolution - its splitting 
into different sub-peaks (Fig. |SJ is in the following ex- 
plained in the framework of the double-resonant Raman 
models The model explains the D' line in the following 
way (see Fig.^a)): An electron is vertically excited from 
point A in the it band to point B in the it* band by ab- 
sorbing a photon. The excited electron is inelastically 
scattered to point C by emission of a phonon with mo- 
mentum q. Since the energy of this phonon (« 150 meV) 
is small compared with the photon energy of 2.33 eV, we 
have drawn the line horizontally, for simplicity. Inelastic 
backscattering to the vicinity of point A by emission of 
another phonon with momentum « q and electron-hole 
recombination lead to emission of a photon with an en- 
ergy about 300 meV less than the energy of the incident 
photon. In principle, two other double-resonant Raman 
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FIG. 5: D' peaks for an increasing number of graphene lay- 
ers along with HOPG as a bulk reference. The dashed lines 
show the Lorentzian peaks used to fit the data, the solid lines 
are the fitted results. The single peak position for the single- 
layer graphene is at 2678.8 ± 1.0 cm" 1 . The peak position of 
the the two inner most peaks for double-layer graphene are 
2683.0 ± 1.5 and 2701.8 ± 1.0 cm" 1 . On the left the value for 
the splitting from double-layer graphene up to HOPG is pre- 
sented. All peaks are normalized in amplitude and vertically 
offset. 



processes, involving the phonons q' and q", are possible 
as well. However, it was argued in Ref. Jlj that their 
weight is very low^ 

In Fig. El we compare the electronic band structure of 
the single layer with the ones of the double layer and of 
bulk graphite. All three dispersion relations were calcu- 
lated from first-principles, using density-functional the- 
ory in the local density approximation, 1 - In the double- 
layer, the 7r and ir* bands split into two bands each. This 
gives rise to four different possible excitations. We have 
calculated the corresponding oscillator strengths 16 and 
found that for the excitation energy of 2.33 eV, tran- 
sitions 1-3 and 2-4 have negligible weight, while tran- 
sitions 1-4 and 2-3 (displayed in Fig.EJb)) have almost 
equal weight. For each of the two dominant vertical tran- 
sitions, there are two possible horizontal transitions. The 
corresponding electron-phonon coupling matrix elements 
for the phonons go to q^ are almost equaltf In theory, 
we therefore expect a splitting of the D' band into four 
peaks of almost equal height. Our experimental data 
(Fig.|5j) shows indeed that the D' line for the double layer 
can be decomposed into four peaks. However, the outer 
two peaks (corresponding to the phonons go & n d gs) have 
very low weight in the experimental data. We calculated 
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FIG. 6: Electronic band structure along the high-symmetry 
lines F-K and K-M: (a) single-layer graphene, (b) double-layer 
graphene, and (c) bulk graphite. For bulk graphite, we display 
the band structure in the direction parallel to the graphene 
planes for different values of the transverse momentum k z . 
Vertical arrows denote vertical transitions by 2.33 eV from 
a valence (n) band to a conductance (n*) band. Horizon- 
tal arrows denote transitions between two states of almost 
equal energy by coupling to a phonon of momentum qi (the 
corresponding phonon frequencies are displayed in Table [IJ • 
Dashed horizontal lines denote transitions with considerably 
less weight than the solid horizontal lines (see text). 



from first principles 1 - the phonon- frequencies v\ and V2, 
corresponding to the momenta gi and g2- The frequen- 
cies of the highest optical branch are given in Table [I] 
Due to the weak inter-layer coupling the degeneracy of 
this branch is lifted. However, the frequency difference 
remains weak (< 1 cm -1 ) and does not significantly con- 
tribute to the experimentally observed splitting about 19 
cm -1 of the D' line (see Fig.[3J). Table[I]furthermore gives 
the value for 2(v2 — vi). We note that the value obtained 
from our first-principles calculation is only half as large as 
the experimentally observed splitting of about 19 cm . 
This discrepancy is related to the fact, that the double- 
resonant Raman model based on ab-initio calculations 
also predicts a value for the dispersion of the D' line with 
incident laser energy that amounts only to about half 
of the experimentally observed value of 99 cm _1 /cV. 18 
We conclude therefore that the double-resonant Raman 
model can qualitatively explain the fourfold splitting of 
the D' line in the double-layer, but the amount of the 
splitting and the relative heights of the peaks are not 
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TABLE I: Frequencies of the optical phonons involved in the 
double-resonant Raman model. The corresponding phonon 
momenta qi and (?2 are determined from the ab-initio elec- 
tronic band structures of Fig. |S| The splitting of the fre- 
quencies in the double-layer and bulk is due to the (weak) 
inter-layer interaction. 



properly described within this models 

In bulk graphite, the 7r and n* bands split into a con- 
tinuum of bands, i.e., they disperse in the direction k z 
perpendicular to the layer. In Fig. Etc), we display the 
bands for three different values of k z . In the joint-density 
of states, the vertical transitions for k z = have the 
dominant weight and are thus considered in our calcu- 
lations. Since the splitting between the bands is much 
more pronounced than for the double-layer, the value for 
2(^2 — v\) is about a factor of two higher than for the 
double layer. This is in agreement with the experimen- 
tal data, see Fig. El where the splitting increases likewise 
by about a factor of two between the double layer and 
bulk graphite. As in the case of the double layer, there 
are quantitative differences between theory and experi- 
ment for graphite as well: First-principles calculations 
of the oscillator strengths and of electron-phonon cou- 
pling matrix elements predict an almost equal height of 
the peaks whereas the experiment shows that the lower- 
frequency peak has a strongly reduced weight. The peaks 
corresponding to the horizontal transitions qo and q% are 
missing altogether in the experimental spectrum. 

Even though some quantitative differences remain, the 
double-resonant Raman model explains well the observed 



differences in the D' line as we go from the single-layer via 
few-layer systems to the bulk limit. The quantitative dif- 
ferences may be an indicator that some essential effects 
are not properly included in the model. E.g., the role 
of quasi-particle effects (electron-electron interaction)^ 
and of excitonic effects (electron-hole interaction) in the 
double-resonance process remains to be understood. The 
importance of these effects has been recently demon- 
strated for electronic excitations in carbon nanotubes 
(both semiconducting and metallic). 21-22 A similar im- 
portance may be therefore expected for processes that 
involve electronic excitations in graphite. 23 

In conclusion, Raman mapping reveals to be a pow- 
erful tool to investigate single- and few-layer graphene 
flakes. It turns out that the width of the D' line is highly 
sensitive to the crossover from single- to double-layer 
graphene, which is explained by a peak splitting follow- 
ing the double-resonant Raman model together with ab 
initio electronic band structure calculations. A remain- 
ing open question is the decrease of the G line intensity 
with decreasing layer number compared to the almost 
constant spectral weight of the D' line and the accompa- 
nied upshift of its frequency for double- and single-layer 
graphene. The structural quality of the flakes is studied 
by analyzing the D line intensity: no defects are detected 
in the inner part of the flake. The D line signal from the 
boundaries of the individual sections of the flake suggest 
that they act as elastic scatterer. 
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